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Nonlinear  Theory  of  a Bearing  Surface  of  Arbitrary  Extent 
by  A.  N.  Panchenkov 

(Institute  of  Hydromechanics,  Academy  of  Sciences,  UkSSR) 

The  nonlinear  problem  of  a wing  of  small  extent  in  liquid 
and  gas  is  one  of  the  difficult  problems  of  fluid  mechanics. 

These  difficulties  are  primarily  related  to  the  complexity  of 
the  phenomena  arising  in  a fluid  when  the  motion  of  a wing  of 
small  extent  moves  with  a large  angle  of  attack.  The  formation 
of  strong  turbulent  eddies  at  the  leading  and  lateral  wing  edges, 
the  large  stream  angle  of  taper,  the  formation  of  a turbulent 
region  beyong  the  wing — all  result  in  the  hydromechanical  character- 
istics of  the  wing  being  nonlinear  with  respect  to  the  angle  of 
attack  so  that  the  nonlinearity,  with  its  attendant  increase  in 
the  lift  on  the  wing,  is  substantial  even  at  small  angles  of 
attack.  This  severely  handicaps  the  mathematical  analysis  of  the 
phenomenonjslnce  the  fundamental  set  of  assumptions  in  the  linear 
theory  of  a bearing  surface  are  not  even  applicable  here. 

In  conjunction  with  the  mathematical  difficulties,  the 
well-known  results  from  the  nonlinear  theory  of  a wing  of  small 
extent  are  obtained  on  the  basis  of  certain  physical  pictures  of 
the  flow,  based  on  experipefrital  facts.  Many  investigators  have 
assumed  various  flow  models  and  have  obtained,  within  the  frame- 
work of  their  models,  results  in  agreement  with  experiment 
(2,  3*  8*  9»  111  • These  results  are  definitely  of  interest 
from  the  point  of  view  of  a basic  explanation  of  the  phenomena 


involved,  the  pattern  of  these  phenomena,  and  many  technical 
applications.  However,  a lag  in  mathematical  investigations  has 
lead  to  a restriction  in  the  class  of  problems  studied  and  to  a 
lack  of  rigor  and  consistency  in  our  concepts. 

The  purpose  of  this  article  is  to  get  theoretical  results 
from  the  nonlinear  theory  of  a bearing  surface  by  mathematical  ana- 
lysis. In  f7]  are  formulated  two  conditions  for  the  "complete" 
approximation  of  a certain  boundary-value  problem  by  means  of  a 
new,  different  boundary-value  problem.  Now  it  is  easy  to  show 
that  the  linear  problem  of  a bearing  surface,  derived  from  the 
nonlinear  problem  by  the  classical  perturbation  technique,  does 
not  yield  a uniformly  convdrgent  approximation.  It  is  know  that 
where  a bearing  surface  moves  in  a fluid  a semi- infinite  region 
is  formed  beyong  it,  in  which  the  velocity  potential  and  its 
derivatives  possess  discontinuities.  If  By  is  the  boundary  of 
this  region,  then 


— oo), 


where  p is  some  avergaed  stream  taper  angle  at  a point  at  infinity 
beyond  the  wing. 

In  all  known  linear  theories (see,  for  example, \k\  t\p\  ) 
the  turbulent  region  is  reckoned  from  the  surface  and 


§ = ° 


but  then  for  the  metric  in  R*  space  we  have 

<f(B«;floi)  = 0(P)|x|  (1) 

and  the  condition  in  £7}  for  large  |x|  is  not  fulfilled.  The 
stream  taper  angle  of  the  wing  has  the  following  asymptotic  values* 


— 0- 


(2) 


Wt- 


p==Q(x) 

P = Q (o)  .(X  - o), 

where  “ is  the  wing  angle  of  attack  and  X is  the  relative  span 
of  the  wing. 

In  this  problem,  it  is  sufficient  to  consider  the  region  £2*. 
jcGI  + co-hCj  where  C is  some  constant  and  then,  for  a<l;.  and 
the  classical  linear  theories  give  a uniformly  convergent  approximation. 

It  must  be  noted  that  perturbation  theory  gives  an  "incomplete" 
approximation.  In  the  well-known  linear  theories,  the  operator 
P = 1 and  the  operator  T = T„.  The  "complete"  approximation  has 
a supplementary  freedom  of  choice  in  the  operators  P and  T0  and 
can  yield  a better  solution.  In  this  sense,  the  method  of 
Poincare-Lighthill-Go  fiolis  quite  interesting  since  it  is  based 
on  the  idea  of  a "complete"  approximation. 

Let  us  consider  the  motion  of  a bearing  surface  of  arbitrary 
extent  with  a constant  velocity  v in  an  unbounded  fluid.  Ljiet 


The  velocity  potential  9 cannot  be  a continuous  function 
throughout  the  region  outside  s,  since  there  will  exist  beyond  the 
bearing  wing  a semi-infinite  region  i^,  in  which  the  velocity 
potential  and  some  of  its  derivatives  will  undergo  discontinuities. 
Indeed,  the  presence  of  a semi-infinite  region  determines  all 
tne  difficulties  in  this  problem.  In  the  region  , the  velocity 
potential  will  even  satisfy  the  Poisson  equation  (5)  but  with  an 
unknown  density  p.  To  get  a solution  with  the  required  properties, 
we  need  supplementary  conditions  determining  the  flow  structure  in 
the  region  As  one  such  condition,  we  can  take  the  physical 

condition  of  pressure  continuity  in  the  transition  of  any  surface 
to  As  examples  of  such  surfaces  we  can  take  surfaces  s^,  on 

which  the  density  p is  constant. 

If  <pSj  is  a tangential  derivative  of  the  potential,  then  the 
condition  of  pressure  continuity  during  transition  of  surface 
can  be  written  in  the  form 

9S+  — 9 - = 0;  Q € s, . ( 6 ) 

Moreover,  the  solution  of  the  problem  must  satisfy  the  supplementary 
conaitions  related  to  satisfying  the  Zhukovskiy-Chaplygin  postulate 
wnich  can  be  written  as: 

9s+~9  =0;  QCL,  (7) 

where  L is  the  equation  for  the  trailing  edge  of  the  wing. 

The  difficulties  of  the  problem  formulated  are  clear.  The 
region  U-^  is  not  known  beforehand  and  also  unknown  is  the  density 
in  the  Poisson  equation  (5).  The  region  in  dimensionless 

coordinates  such  that  the  wing  span  is  equal  to  1,  has  a thickness 

& = o(y)  and  thus  we  obtain  one  of  the  required  conditions  for  the 

linearization  of  the  problem: 


T « 1 


(8) 


Under  condition  (8),  the  region  can  be  approximately  replaced 

by  a certain  surface  s^.  Well  known  linear  theories  [4,  5]  start 
out  from  the  assumption  that  the  surface  s^  is  inclined  only  slightly 
from  the  xOy  plane  and  condition  (7),  having  the  form: 

<PX+  — <PX-  = 0, 

is  satisfied  on  the  plane  s , lying  in  the  xOy  plane. 

p A 

We  shall  solve  the  problem  posed  by  using  the  idea  of  a 
"complete"  approximation.  Let  us  introduce  a new  three-dimensional 
space  R^  with  coordinates  x1,  y^,  z^.  We  take  the  relation 

3 3 

between  the  coordinates  in  the  and  R^  space  in  the  following 
form: 

x = *i  + Fx  (at,;  yx\  z,); 

y=  yi  + Ft  (x,;  yx,  2i);  ^ J 

z^Zi+F,  (*,;  yx\  Zj), 

where  the  functions  Fj  belong  to  the  class  FfeC^Q).  Let  the  norms 
of  the  functions  F.  and  their  first  and  second  derivatives  in  the 

2 J 

class  C have  the  values: 


\\F,\U  = 0(t).  (10) 

where  e is  some  small  parameter  in  the  problem.  Then,  when  the 
conditions  (10)  are  fulfilled,  the  potentials  <p  and  space 
coordinates  can  be  written  as  a series: 


9 = <P«  (•**•  yi>  2i)  + eq>,  (x„  yx,zj  + e2<ps  (x„  yx,  z,)  + 

X/  x*j  ex)  (x,,  yt,  Zj)  + etxj(xl,  y, , z(); 
x,  = x;  yt  = y;  x,  = z. 

3 

In  the  zero ' th  approximation  in  the  R^  space,  the  potential  s0 
will  also  satisfy  the  Laplace  equation.  In  the  usual  formalism, 
the  problems  for  the  functions  <p/(£,  t),  £),  can  also  be  treated, 
wherein  they  will  even  satisfy  the  Poisson  equation.  Since  it  is 
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clear  from  the  physical  nature  of  the  problem  considered  that  the 
greatest  influence  on  the  wing  characteristic  will  be  exerted  by 
the  presence  of  the  stream  taper  angle  at  infinity  beyond  the  wing, 
we  shall  limit  ourselves  in  what  follows  to  the  problem  involving 
deformation  of  only  the  y coordinate.  Thus,  we  shall  set: 


x — 

y=  y,  + ey\(xvyr  z>)  + zi): 

2 = Z,. 


(11) 


In  this  method,  there  is  a freedom  of  choice  of  the  functions 
<P(/» (I, Tl. £).  which  should  be  so  chosen  that  they  meet  the  approximation 
conditions.  According  to  the  second  approximation  condition  (b) 
in  reference  [7]: 


£ = 2 + j <Pcdx 


(12) 


and  to  a first  approximation: 


£ = 2 + ] <Poc dX 


(13) 


Let  us  introduce  into  consideration  the  function  Q=— For  the 
function  0 in  the  space  we  have  the  following  problem: 


A©  = 0;  ‘ 

e„i  = Pi  (giY  gi  € Sul 

0-0;  x -*  + oo; 
e+  — e_  = 0 na  L,; 

0 = 0;  g,  € SXI. 


(14) 


Let  us  take  the  solution  of  the  problem  to  be  in  the  form: 

Q (P)  = $ Y (g)  K (p.  g)  dg<  (15) 
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where  K(p,  g)  is  the  fundamental  solution  of  problem  (14)  correspon- 
ding to  a dipole. 


In  the  linear  theory,  in  (15)  they  take  instead  of  s^. 
In  this  theory  we  are  developing,  it  is  sufficient  to  limit 
ourselves  to  the  approximation 

Q (p)  ~ jy  (gi)  Ki  dgt\ 

hp 

K (Pigi)  = K l(*i  — Si);  (yi  — ^i);  & — £));  ( 1 6 ) 

Ki  (Pigi)  = K[(x  — ( y — tij  (z  — SO  + F (Si)] . 

where  s^  is  the  projection  of  the  bearing  surface  onto  the  x^Qy^ 
plane.  Continuing  on,  we  have 

Qm  ~ Qi  — Pi  (gi)'<  gi  € s„; 

51 P 


(17) 

(18) 


Going  over  to  the  potential  % and  carrying  out  the  computations, 
we  get  from  condition  (17)  the  following  integral  equation  for 
the  problem: 


1 ° , y — r) 

<P#*  = 4n  J m dyAto  — t))2  + F2  (|)  ' 

s1p 

\(yl-r]f  + {xl-lY  + 21^(1)]  1) 

(*-S)K(y-T»2-M*-&)2  + /'2(£)  J|  gl 


(19) 


For  further  work  it  is  convenient  to  switch  to  dimensionless 
coordinates 


y_ . 
B • 


M*0  = 


B 

c<nr 
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where  C(n)  is  the  wing's  semi-chord  in  cross  section  r\  and  B 
is  the  wing's  semi-span.  Then 


1 [X2  to)  (y  - T -f  (x  - £)*  + 2P  (l)  I 

(x~l)  + (*-%)* + F(i) 


Equation  (20)  solves  the  problem  of  the  motion  of  the  bearing 
surface  of  any  arbitrary  shape  for  average  attack  angles  a. 

The  function  F(|)  is  determined  by  the  nature  of  the  transformation 
in  (12)  in  the  region  of  the  bearing  surface.  For  F(l)  =0,  equation 
(20  yields  well  known  results.  The  function  F{1)  from  the  point 
of  view  of  the  flow  mechanics  determines  the  distance  of  the 
vortices  from  the  bearing  surface  at  any  point  . 

Since  the  nonlinear  effects  are  most  substantial  on  wings  of 
small  extent,  we  shall  study  more  closely  the  problem  of  a wing 
of  small  extent.  For  A->-0,  equation  (20)  takes  on  the  following 
correct  limiting  form: 


Let  us  consider  a solution  for  the  function  F(Q: 

f(S)  = (fc-l)PlT))  (22) 

Fulfilling  condition  (4)  on  the  trailing  edge  of  the  wing,  we 
have,  for  the  function  ori  > the  equation 
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r 


I 

2 n\ 


+i 


(23) 


C %(y  — ' i) 


sin  a. 


Function  m-  in  the  form  (22)  and  equation  (23)  encompass  a 
large  group  of  nonlinear  physical  theories.  In  connection  with 
this,  we  note  one  important  fact  related  to  the  choice  of  angle 
g.  It  is  easy  to  show  that  for  the  angle  g,  the  correct  asymptotic 
value  is: 


P = 0 (Xx)  (x  > 0).  (24) 


The  value  in  (24)  is  obtained  from  equation  (23)  under  the  condi- 
tion that  the  norm  of  a>  is  bounded  as  A+0 . If  || <D ||x-o=  M (where  M 
is  some  number),  then  x=0.5*  In  all  known  theories  [2,  3,  5,  7, 
10],  the  angle  g is  taken  as  finite.  In  a series  of  works  on  this 
problem,  it  is  taken  as  equal  to  half  the  flow  taper  angle  and; 

in  other  works,  various  other  means  are  employed  for  its  computa- 
tion. Of  course,  it  is  basic  that  for  finite  values  of  the  angle 

g (A-»-0),  the  results  of  the  theory  must  give  infinite  values  for 

the  lift  on  the  wing  — a prediction  which  does  not  agree  with 
experimentally  observed  results. 


Equation  (21)  makes  it  possible  to  obtain  detailed  information 
on  wings  of  small  extent  of  arbitrary  shape.  The  determination  of 
the  function  F(l)  is  quite  important.  Since  it  is  not  evidently 
possible  to  completely  describe,  within  the  framework  of  the  theory 
of  an  ideal  fluid,  the  physical  phenomena  during  motion  of  a wing 
of  small  extent,  the  function  F(£)  makes  possible  a more  precise 
reckoning  of  the  influence  of  viscosity,  if  in  determining  this 
function  certain  supplemental  hypothesis  are  introduced  concerning 
the  flow  pattern.  However,  in  this  regard  it  must  be  noted  that 
even  without  the  attraction  of  supplementary  data  on  the  influence 
of  viscosity  within  the  framework  of  an  ideal  fluid  results  are 
obtained  agreeing  satisfactorily  with  experimental  data.  To  get 
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further  results,  let  us  take  the  angle  3 in  the  form 


P = a-ai,  (25) 

where  on  is  the  flow  taper  angle  at  infinity  beyond  the  v/ing. 

In  connection  with  definition  (25)  and  the  evaluation  in  (24), 
still  one  more  important  assumption  must  be  noted:  formulas  (24) 
and  (25)  assumed  the  presence  of  more  information  than  that  given 
by  the  linear  theory  for  a wing  of  small  extent.  The  basic  result 
of  the  theory  of  a wing  of  small  extent  (Namely  that  the  flow 
taper  angle  equals  the  angle  of  attack)  leads  to  the  result  that 
3=0  whereas  we  need  information  of  the  type  in  (24). 

If  we  make  the  Laidlow  approximation  for  a wing  of  arbitrary 
extent  (11)  and  determine  the  flow  taper  angle  from  the  theory 
with  arbitrary  extension,  then  we  get 


XC.(X)  + 2[C,(X)-  1], 
P XC2  (i)  + 2 [C,  (1)  +1]  * 


(26) 


where  C^(X)  can  be  determined  from  the  work  of  Nishiyama; 


Cl(*)=  1 -< 


/X2  + 4 


(27) 


for  A->-0,  we  have: 


2]/2 


(28) 


The  results  in  (28)  coincides  with  the  value  in  (24)  for  x = 0 . 5 . 
Equations  (21)  and  (23)  can  be  solved  by  known  techniques  [5]. 
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To  get  finite  results  for  the  lift  coefficient,  the  solution 
to  equation  (23)  is  obtained  via  a variational  approximation.  The 
lift  coefficient  is  obtained  in  the  form: 


C„  = 


jtX 


' 2 ’ 


(29) 


C - — * / 

^6  = 2 je  * - 


lUk) 


dk. 


The  function  w is  calculated  in  the  form: 


= 0,  5t|  + 0.  25x4  + 0,  0625-r®  + 0,  0469r«  + 
+ 0.  0237tg°  + 0,  0188t'2  + 0,  098  It'4; 


*+i_t 

^ 2X' 


(30) 


In  conclusion,  let  us  consider  the  problem  of  a wing  of 
arbitrary  extent.  In  this  problem,  we  shall  introduce  into 
equation  (20)  the  approximation: 

l*-2  (y)  (y  - -n)2  + (*  -£)2  + 2P-(S)]  _ 


Vw  (y)  (y-  ti)2  + (x  - if  + F2  © 


; Cj  (X)  \x  — II  + c2  (X)  X (T))  sign  If/ — T||- 


(3D 


Then  equation  (20)  goes  over  into  the  following  form: 


<P*.  = 


" 2n  J 


F2(x) 


+i 


(32) 


4;t 


11 


F 


I 


P 


The  approximation  equation  (32)  is  significantly  simpler 
than  equation  (20)  and  can  be  effectively  solved.  The  following 
step  is  a transition  from  equation  (32)  to  an  equation  of  the 
Prandtl  type: 


r(ir) 


2 X(j/)C4(X) 


a(y) 


C,(M 
2 n 


T'  (ti)  (y  — ti) 


(y  — »i)2  + 


fJl 

(X  ft)) 


(33) 


where 


C3  (X)  — 1 + Cj  (X); 

F = F(—  1);  r(T))  = -i-’jvte.)  dl 

— I 

For  F=0  and  CgtX^C^A)*!,  equation  (33)  becomes  the  well-known 
Praddtl  equation  for  a wing  of  large  relative  span. 

The  functions  C^(A)  and  C^(X)  should  be  selected  from  the 
condition  of  best  approximation  in  (31)*  Equation  (33)  makes 
it  possible  to  study  a wide  class  of  problems  posed  by  the 
Prandtl  equation  in  linear  theory. 


As  a first  result  we  obtained  formulas  for  the  lift  force  and 
parasitic  drag  of  a wing  with  an  elliptic  circulation  distribution: 


C„  = 


C2  (X)  1 + 


‘-AW 

C2  (X)nX 


V 


cr  = 


C 

nX  e 


(34) 

(35) 


As  follows  from  (35),  the  wing's  parasitic  drag,  obtained 
from  the  nonlinear  theory,  has  turned  out  to  be  less  than  that 
obtained  from  the  linear  theory.  Formulas  (34)  and  (35)  are  exact 
in  the  linear  theory,  but  in  the  theory  being  considered  they 
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are  obtained  via  a variational  approximation,  since  an  elliptical 
distribution  does  not  even  allow  an  exact  solution  of  equation  (33). 


1 

As  is  known  from  experiments  [11],  for  large  attack  angles, 
some  flow  regimes  for  a wing  of  small  extension  are  possible 
depending  on  the  shape  of  the  j.eadjng  and  lateral  edges  of 
the  wing.  For  a wing  with  rounded  edges,  there  is  observed  a more 
uniform  distribution  of  the  rising  vortices  along  the  wing 
surfaces,  while  for  sharp  edges  there  is  observed  strong  concen- 
trated vortices  at  the  lateral  edge,  with  tne  result  that  the 
lift  in  the  latter  case  is  significantly  higher  than  that  for  a 
wing  with  rounded  edges.  Function  F(£)  from  (22)  and  angle  0 
from  (26)  and  (28)  correspond  to  a wing  with  rounded  edges. 


Fore  the  second  regime,  one  must  take  another  distribution  for 
the  function  F(£,  n)  along  the  wing  surface.  From  experiment  it 
is  known  [11]  that  the  lift  with  sharp  edges  is  approximately  two 
times  larger  than  that  for  a wing  with  rounded  edges.  To  a first 
approximation,  we  shall  characterize  this  regime  by  the  angle: 


(36) 


Calculation  of  the  functions  V.  and  C 

p y 

(30),  and  (3^)  with  angle  0 from  (36)  gives 
with  known  experimental  data  [8,  11]. 


from  formulas  (29), 
good  correspondence 
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